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We present a first theoretical study of a merging of diabolical points in the context of supercon- 
ducting circuits. We begin by studying an analytically solvable four-level model which may serve 
as theoretical pattern for such a phenomenon. Then, we apply it to a circuit named Cooper pairs 
pump, whose diabolical points are already known. 
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I. INTRODUCTION 

Superconducting circuits^ become more and more cen- 
tral in modern quantum physics. Their principal build- 
ing blocks are ultra-small Josephson junctions^ which 
can be assembled in a variety of different ways, each 
of them defining a specific quantum system. Often eas- 
ier to manipulate, these solid state devices tend to sup- 
plant atomic and photonic systems, the "old paragons" 
of quantum mechanics. Notably, superconducting cir- 
cuits are widely used to engineer qubits'^^^: the non- 
linear behavior of Josephson junctions serves to isolate 
couples of levels in a Hamiltonian spectrum. They are 
also used to perform analogs of cavity quantum electro- 
dynamics^"^^ (a qubit plays the role of an artificial atom 
while a transmission line carries artificial photon modes) , 
(non-) Abelian holonomies"'^^""'^^, (non-) Abelian quantum 
charge pumpings^*""^^, etc. In brief, they are good candi- 
dates to implement quantum logic operations'"'^^ and also 
promising for applications in electrical metrology^^. The 
Cooper pairs pump (CPP) which we will consider in this 
article is an archetype of quantum circuit having a few 
(collective) degrees of freedom. It has been demonstrated 
theoretically in Ref. 22 a possible topological quantiza- 
tion of the pumped charge, through an invariant called 
first Chern number (or Chern indcx)^^'^^. It relies on the 
existence of diabolical points^^ in the three-dimensional 
parameter space of the system, i.e. on double degenera- 
cies characterized by a linear dispersion in all directions 
of that space. 

Relatively recently, G. Montambaux et al. have 
demonstrated the possibility of merging Dirac points in 
certain two-dimensional crystals, especially in hexago- 
nal — graphene-like — lattices^^"'^^ (see also Refs. 33 
and 34). Dirac points are nothing else but diabolical 
points in the reciprocal space of crystals. They are "nat- 
urally" located at points of high symmetry — e.g. at 
vertices of a regular hexagonal lattice. But, in accor- 
dance with the famous Wigncr-von Neumann theorem^^, 
they may move, driven by well-chosen additional param- 
eters. In the graphene example, the two triangular sub- 
lattices carry non equivalent Dirac points. The merging 
of two neighboring non equivalent Dirac points evokes 
the meeting of a knot and its anti-knot: monitored by 
a merging parameter, they move closer together, then 
merge into a single degeneracy and finally disappear. At 



the transition, the single degeneracy is characterized by 
a quadratic dispersion in the merging direction. Inspired 
by the works of G. Montambaux et ai, we present in this 
paper a theoretical study of a merging of diabolical points 
in the context of superconducting circuits. The choice 
of the CPP was motivated by its well-known diabolical 
points located in a hexagonal lattice, a property which 
confers to that circuit a great similarity with graphene. 

This paper is organized as follows. Section II is devoted 
to the short study of a merging of diabolical points using 
a four-level model. It is interesting by itself since it pro- 
vides a "universal Hamiltonian" which is susceptible to 
be realized in diflFerent quantum contexts as ours. It gives 
en passant the opportunity to briefly review some char- 
acteristics of double degeneracies in a parameter space. 
In Sec. Ill, we present the CPP and emphasize the sym- 
metry origin of its "mobile diabolical points". Finally, 
via a modification of the circuit, we suggest in Sec. IV a 
way of merging these points. This will be done through 
an effective Josephson energy as merging parameter. 



II. THE FOUR-LEVEL MODEL 

We consider a model Hamiltonian depending on a 
triple of tunable parameters R = {X, Y, Z) and having 
the form 
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in an orthonormal basis {|ei), |e2), jea), \e4)}. Here, F ^ 
and ^ are constants, and G is an additional tunable 
parameter. The latter is dubbed merging parameter for 
a reason which will appear shortly. The set of vectors 
H forms the natural parameter space of the problem. 
In this space, the spectrum of H possess the symmetry 
T>2h- Indeed, H is (anti)unitary transformed under sign- 
reversing oi X ,Y and Z . Explicitly, we have 

1. H{X,Y,-Z) = JCH{X,Y,Z)J0 where JC is the 
complex conjugation operator with respect to the 
basis {|ei), |e2), \e^), \ei)}\ 



(a) " (b) (c) " (d) 

FIG. 1. Plots of the two lowest energy levels as functions of X and Y, with Z = and ^ = for different values of the 
merging parameter G. As long as G > Gc the distance between the diabolical points and D~ decreases with G. They 
merge for G = Gc into a single degeneracy characterized by a quadratic dispersion 5{X) in the X-direction and disappear as 
G < Gc. The plot unit is |F|. (a) G = Gc + 0.25 \F\; (b) G = Gc + 0.05 \F\- (c) G = Gc; (d) G = Gc - 0.3 \F\. 



2. H{X, -Y, Z) = [U{Z)K] H{X, Y, Z) \U{Z)K\ ^ with 
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3. H{-X,Y,Z) = [U{Z)T]H{X,Y,Z)\U{Z)T]'^ with 
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In particular, under ii-inversion, wc observe the very 
simple unitary equivalence H{—R) = TH{R)T^ . 

Since (ei|(iJ - fi){H - i^)\e4} = 2f 2g2^z ^ g, for any 
fi and v, H has at least three distinct eigenvalues. Thus, 
A is a (doubly) degenerate eigenvalue of H if and only if 
(iff) it exists a real (3 > such that 



{H - X){H + A + I3){H + A - /?) = 0. 



(2) 



In this case, —A — /3 and —A + /3 are the other eigen- 
values and A is the smallest one iff /3 < — 2A. The ma- 
trix element of Eq. (2) between (ei| and |e4) imposes 
A = — ^ — GcosZ. This value is a root of iJ's char- 
acteristic polynomial iff one of these three conditions is 
verified: (i) F = G = ; (ii) X = ±{2^ + GcosZ) ; (iii) 
y = and Z = mod tt. A study of Eq. (2) within 
conditions (i) or (ii) shows that they can not give rise to 
a ground degeneracy. 

We now consider the last possibility. It is sufficient to 
restrict ourself to the case Y ~ Z = since Z ^ Z + tt 
amounts to change the sign of F and G. After a little 
algebra, we find that the ground level of H is degenerate 
iff G is greater than the critical value Gc = y^2F^ + £^—£^ 
and X = ztXd with 
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(G + 2e)(G + Gc + 2e). 



As long as G > Gc, the two distinct points D^, lo- 
cated at = (=F^d, 0, 0), are isolated degenerate points 
in the i2-space. They illustrate a classic theorem of von 
Neumann and Wigner'^^ which states that, generically, 
twofold degeneracies have codimension three. Alterna- 
tively stated, in an A'^-dimensional parameter space, it 
exists generically submanifolds of dimension iV — 3 over 
which a level is doubly degenerate. Here, = 3 and the 
degenerate subspaces are points (dimension: 3 — 3 = 0). 
One should also think about a bigger space collecting 
all the parameters entering in H, as the six-dimensional 
space of vectors ^ = (A, y, Z, ^, F, G). The subman- 
ifold verifying simultaneously F ^ 0, Y = Z = 0, 
G > Gc(^,F), and X = ±Xd{^,F,G) carries a degen- 
erate ground level and has a (co)dimension 3 in the 
space, as expected. 

The points move closer together when we reduce 
G, while conserving the symmetry relations [H,U] — 
[H,IC] = (see Fig. 1). It corresponds to a generic sit- 
uation: if they would deviate from the planes y = or 
Z = 0, each of them would split into 2 (or 4) distinct 
degenerate points. They merge at ii = for G = Gc 
and finally disappear as soon as G < Gc, the minimal 
gap between the two lowest levels being 



A - - i^^l6F^ + {G~2^y^^3G 



Let n be the projector into iJ(i?+)'s ground eigenspace 
and {|1), 1 2)} an orthonormal basis of that subspace. Let 
c = {o'x ,<^y,<^z) be the triple of operators represented by 
the usual Pauli matrices in the basis. Up to an unim- 
portant component along II, it exists a unique fixed real 
nratrix M of order 3 such that n\VH{R+) ■ {R- R+)]ll 
reads er ■ M{R — il+). Obviously, M depends on the 
choice of the basis, but it is a simple task to show that 
the signum of its determinant is intrinsic to the degener- 
acy. It is the signature^^ of the degenerate point in 
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the i2-space. One can explicitly choose 
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Within this choice, one finds 
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As long as G > Gc, the signature of L>+ is +1. In 
its vicinity, the two lowest levels are close together and 
(7 ■ M{R— i?"*") is an accurate Hamiltonian for the states 
belonging to them. Since has a nonzero signature, 
the level splitting around it is effective from the first or- 
der in JS+II in all directions of the parameter space: 
is a diabolical point^^. By symmetry, so does D~ , 
whose signature is found to be —1. For G = Gc, the 
single degeneracy located at the origin has a vanishing 
signature, because the dispersion in the merging direc- 
tion is quadratic. A perturbativc analysis shows that a 
deviation (0, 0, 0) — !> {X, 0, 0) opens a gap 
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Before introducing the system which will serve to re- 
alize our four-level model, let us end this section by two 
remarks. The first one is peculiar to the model: if y, Z 
are suppressed and ^ F tunable, the Hamiltonian may be 
used to construct non-Abelian holonomies"^^'"^^ over the 
manifold satisfying simultaneously F 7^ 0, G Gc{^,F) 
and X = ±Xd{^, F,G). The model may also serve to 
implement non-Abelian pumpings having Z as pumping 
parameter^^. The second remark is more general and 
concerns the signature. Consider some Hamiltonian H 
continuously defined over a ii-space. Suppose the exis- 
tence of a nonsingular transformation t : R t-^ R' as- 
sociated with a fixed symmetry operator T, such that 
H{R') ^ TH{R)T'^. If T is unitary, it is straight- 
forward to verify that the signatures are conserved by 
the transformation if t is orientation-preserving and re- 
versed otherwise. If T is antiunitary, the contrary oc- 
curs. In our example, the signature of D~ is due to 
the orientation-reversing map R ^ —R associated with 
the unitary operator T. Moreover, we will below use 
successive orientation-preserving transformations of the 
parameter space, without incidence on the signature. 




FIG. 2. The Cooper pairs pump (CPP) is essentially an array 
of three Josephson junctions in a loop configuration. It de- 
pends on three external parameters: the gate voltages V^k on 
the superconducting islands and the magnetic flux t&x thread- 
ing the loop. The system is said mirror symmetric if the "ex- 
terior junctions" are identical (i.e. i?jL = -Ejr, Cl = Cr) and 
totally symmetric if all the junctions are identical. 



III. THE COOPER PAIRS PUMP AND ITS 
DIABOLICAL POINTS 

A. Basic settings 

We consider one of the simplest implementation for a 
CPP, represented in Fig. 2. It is a small- inductance su- 
perconducting loop (F — ^ 0), threaded by a magnetic fiux 
$x and broken by three ultra-small Josephson junctions2. 
The junctions arc assumed non-dissipative. They enclose 
two superconducting islands, polarized by gate voltages 
Vgk through low gate capacitances Cgk {k ~ 1,2). We set 



'-'i^'-^i^'' the corresponding gate charges in unit of 
2e (> 0). The vanishing loop inductance leads to a bias- 
ing phase2 = 27r-|^ accross the CPP, where '^o = ^ 
is the quantum of flux. We suppose $xj l^i and Vg^ 
independently tunable. The system has two collective 
degrees of freedom — one for each island — and de- 
pends on three external parameters that we choose to 
be rigi, ng2 and if-x_ rather than V^i, and $x- The 
conjugated operators assigned to the degrees of freedom 
are and ipk [k = 1,2): the number of Cooper pairs 
in excess (with respect to charge neutrality) and the 
phase of the superconducting parameter of the fc"^ is- 
land, respectively. They verify the commutation relations 
[rij, (y9fe] ~ i Sj^k- Wc will study the CPP in charge repre- 
sentation; \ni,n2) will represent the fundamental charge 
states defined such that e^"^^|ni,n2) = \ni =F l,n2) and 
e^"^^|ni,n2) = |ni,n2 T 1)- Since eigenvalues of ni and 
71,2 can theoretically cover all the range of Z, the orthonor- 
mal basis ^ = {|ni, 712) | [111,112) G ^2} spans the whole 
Hilbert space of the problem. 

Throughout this article, Ge = Gl + G + Gr will be 

the capacitance unit and Eq = -^^^ the energy unit. 
The latter is a typical charging energy of the circuit. We 
will study the system in the Coulomb blockade regime. 
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characterized by Josephson energies small in comparison 
to Eq = 1. Using the canonical quantization procedure, 
a Hamiltonian H = H{ng-^,ng2, (fix) may be derived for 
the system. It splits into two parts: a charging Hamil- 
tonian He ~ Hc{ngi,ng2) and a Josephson tunneling 
Hamiltonian Hj = Hj{ipx). Neglecting the gate capaci- 
tances in comparison to Cs = 1 and using the notation 
a = (ai, a2), the former is 

Fc = ^(n-ng)-C~^(n-ng), (3) 
where C is the capacitance matrix: 

^_ f c\ + c -c \ 

\ -C C + CnJ- 

The charging Hamiltonian is obviously diagonal in the 
basis and verifies i?c(ng + a) = e''"- '^ Hc{ng) e'"- '^ 
for any integer vector a. Over the rig-planc the en- 
ergy surface of the eigenstate |0) — |0,0) is an elliptic 
paraboloid centered at rig =0. Thus, the energy surface 
of \n) = |ni,n2) = e^*" '^|0) is simply the translation 
by n of this paraboloid. Two different states \n) and 
\n') are degenerate on a straight line characterized by 
n, n' and the capacitances. Then, one easily checks that 
\ni, n2) is the ground state of Hq in a hexagon hex(ni, 712) 
centered at rig = n. That defines the well-known hon- 
eycomb lattice of the GPP. It is graphically obtained by 
integer translations of two noncquivalent lattice points 
whose coordinates are 

^ 4 c ( (^:;);; ) ^ 

The lattice picture is useful if we identify each fun- 
damental state |ni,ri2) with its corresponding hexagon 
hex(ni,n2). With respect to _ffc, the common side of 
two neighboring hexagons is a piece of the degeneracy 
line between the states, while the vertices are points of 
triple degeneracy. Introducing the distance induced by 
the scalar product {x\y) ~ 2~^l^x ■ C^^y in the plane, 
this picture allows to interpret the charging energy of 
\ni, 712) as the squared distance between rig and the cen- 
ter of hex(ni, 71,2). "Branching" i/j, which can be brought 
into the form 

H,^ = C/((/3x)| - -EjL cos(<pi + ipx) ~ E,m cos(v?2 + ^x) 

-E]Cos{ip2-(pi+^x)]^U{(px)\ (4) 

couples the neighboring states and lifts generically the 
degeneracies oi He- Explicitly, U{(px) — e"^ '""^", with 
ni = 1- Cr(C-1)i2 and K2 = 2 - Cr(C-1)22. In the 
Coulomb blockade regime, Hj is seen as a perturbation 
of He- As a good approximation, the Hilbert space may 
be reduced to its subspace spanned by a few number of 
fundamental states in the neighborhood of rig. To this 
end, we only take into account the states \ni,n2) at a 
distance of n„ shorter than a certain value'^^ . 



ng2 




FIG. 3. The honeycomb lattice of the GPP over the rig-pIane. 
It is generated by the triple points whose coordinates are 
given in the text. Each state |ni, 722) is the ground eigenstate 
of He in the hexagon hex(7ii, 712). In the space of vectors 
r = {rig, ifx), the spectrum of the Hamiltonian H possess the 
symmetry €21-1 ■ It becomes = C2h x o"^ if the CPP is mirror 
symmetric (the vertical reflection planes are represented by 
thick lines). 



Since i?(ng + a,(px) = e''"- '^ H{ng, Lp^) e'"- '^ , trans- 
lations of lattice vectors a leaves the physics unchanged 
up to a displacement |n) — > \n + a) of the fundamental 
charge states. Moreover, performing the gauge trans- 
formation |rti,rt2) U{ipx)\ni,n2), the Hamiltonian is 
invariant under the translations V'x ^ V'x + 2A:7r (k £ Z). 
Thus, the spectrum of H possesses the translational sym- 
metry of a hexagonal prism lattice in the space of vectors 
r = (rig, ifx). In the new representation, let us introduce 
the complex conjugation operator /C, the "sign change 
operator" S : |7ii,7i2) n- | —711,-712) and the "charge 
exchange operator" V : |7ii,7i2) 1— > \n2,ni). Taking the 
tpx-a-xis vertical, the spectrum of H possesses the point 
symmetry C2/1 = (Jh x i'- under the reflection and the 
inversion l, we have H{r) = lO H{ahr)K. = S'^H{Lr)S. 
In particular, the symmetry C2 = t o Uh implies that 
the Hamiltonians at {T^,ipx) and {T~,ipx) are antiuni- 
tary equivalents, and even unitary equivalents iS ip^ ~ 
mod TT. 

If the CPP is mirror symmetric, that is to say, 
if the "exterior junctions" are identical (i?jL = ^'jr 
and Cl = Cr), the symmetry 2?2ft = C2h x ay is 
reached. Indeed, the reflection ay, shown in Fig. 3, 
exchanges 77,gi and ng2, inducing the transformation 
Hiayr) = [Vi^x)VIC]H{r)[V{ipx)riC]\ where Viip^) = 

g-2j(ni+n2)v5x 



B. The diabolical points 

The CPP is said totally symmetric if the three junc- 
tions are identical. In this specific case, one has 
ng(T^) = ±(i, i). As shown in Fig. 4, the orientation- 
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FIG. 4. The regularized honeycomb lattice of the totally sym- 
metric CPP. The new parameters X and Y are notably chosen 
so that T'^ is located at the origin. 

preserving map (ngi,ng2) {X,Y), such that X = 

^{ng2 + — |) and Y — 5("^g2 ~ ^gi)j makes the 
hexagons regular in tlic {X, y)-plane and places the ori- 
gin at T"*". The charging energy of a state |ni,ri2) be- 
comes the usual squared distance between (X, Y) and the 
center of hex(ni,n2). 

Setting Z ~ (fx, the spectrum possesses the sym- 
metry 2?3/i = a-u X C3 X ah in the so-defined i2-space 
(see Fig. 4). The rotation C3 induces the transfor- 
mation HiCsR) = [V(^x)7^]i^(i^)[V((px)7?.]^ with TZ : 
1 711,712) i-> |1 — 7ii — 772, 77i). The two symmetry oper- 
ators V and TZ generate an unitary representation F of 
the group = C3 x ay in the Hilbert space, such that 
F(cr^) = V and ^{€3) = TZ. At the high symmetry points 
(T"*", (ySx = mod tt), the Hamiltonian is invariant under 
I?3. From the non-abelianity of 2?3 and the faithfulness 
of F, it exists necessarily degeneracies in the spectrum 
of H{T^ ,(px), when ipx = mod tt. More precisely, 
admitting only irreducible representations of degree 1 or 
2, the symmetry allows only double degeneracies (see the 
character table I). 

The ground eigenspace of Hq (T^) — spanned by |0, 0) , 
|1, 0), and |0, 1) — is an invariant subspace of F. If Fg is 
the restriction of F to this space, we have 

/100\ /00l\ 

Fg(o-,„) =001 and Fg(C3) =100 

yoioy yoioy 

in the basis {|0, 0), |1, 0), |0, 1)}. These matrices do not 
commute, so Fg is a non-Abelian representation of ^3. 
Hence, it contains necessarily the two-dimensional irre- 
ducible representation (irrep) E. The trace of Tg^a^) 
shows that Fg contains also the totally symmetric irrep 
A. Obviously, the state 

|yl> = -^(|0,0) + |l,0) + |0,l)) 



TABLE I. The character table of ©a. There are three irre- 
ducible representations: A (totally symmetric), B (antisym- 
metric) and E (two-dimensional). 
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belongs to A. Then, we complete a basis of Fg by choos- 
ing two orthonormal states belonging to E: 




In the basis {|A), \Ei), |i?2)}, we have Tg = A® E with 
A{ay) = A(C3) = (1) and 

By the invariance of i?j(0) and Hj{'k) under any 
Josephson coupling between the states \A), \Ei) and 
|i?2) is forbidden at the high symmetry points. Since 
{Ea\Hj{TT)\E^) - {A\Hj{tt)\A) = < (a = 1,2), up 
to the first order in the ground level of H{T~^,7t) 
belongs to E while the first excited one belongs to A. 
The contrary occurs for H{T'^, 0). Thus, the half-fluxoid 
condition $x = ^ mod <I>o guarantees the double de- 
generacy of the ground level at the point T+. The same 
conclusion holds at T~ from the equivalence between 
H{T-,^^) andi/(r+,(px). 

Let us analyze the signatures of the degenerate points 
= (r±,7r). Redefining, for convenience, .Z as <px — 
is located at the origin of the new i?-space. Using 
the same notations as in Sec. II, basis states |1) and |2) 
of H{D^ys ground level are partners of the irrep E = 
IIFII. They may — and they will — be chosen so that 
the matrices of E are given by (5) in the basis {|1), |2)}. 
Thereby, we have \a) = \Ea) + 0{Ej), a = 1,2. The 
(anti)unitary transformations of H under the action of 
V^fi imply the existence of two reals a and /3 such that 

in the basis {|1), |2)}. The coefficient a may be calculated 
up to the first order in £;j as a = (l|ViJ(0) • 1^12) = 
{Ei\VH{0) ■ -V\E2) + 0{Ej) = 3-1/2 o(£;,j). As for 
13, we find 13 = i(2|Vi?(0) • Z\\) = Z-^l'^Ei 0{El). 
Fixing Z = 0, the two lowest levels, plotted around 
as functions of X and Y , form locally a right circular 
double cone whose aperture is 2arccot(a) = ^ + 0(£'j) 
[see Fig. 5(a)]. This is a consequence of the TD^ symmetry. 
Writing the right-hand side of Eq. (6) in the form cr • 
Mi?, one has dot M = a^/? > 0: the signature of Z)+ 
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(b) 



Cgi 



^jE,C/2 



EI 



^ dc SQUID 

E3-E,C/2 : Cg2 



£;jE,C/2 



Vg2 

iSjE,C/2 



FIG. 5. (a) Plot of the two lowest energy levels as a func- 
tion of X and Y at (^x = tt in the close vicinity of . In 
a totally symmetric configuration, they form locally a right 
circular double cone whose aperture is 2i9 = ^ + 0(i5j). (b) 



Displacements of two diabolical points D and D 
decreases: they merge at their midpoint /. 



as Ei 



is +1. Returning to the natural r-space, the points 
are located at = (ng(r^),7r). Since r = C2r^ and 
H{C2r) = [SlC\H{r+)[SlC\\ the signature of D- is -1. 

The existence of "signed degeneracies" is fundamental 
for the physics of quantum pumpings. They quantize the 
pumped charge along classes of cycles in the parameter 
space^^. They are robust in the sense that their exis- 
tence is guaranteed by the Wigner-von Neumann theo- 
rem even though the T)^ symmetry is broken. In some 
ways, one should say that the symmetry plays the impor- 
tant role of producing signed degeneracies which become 
"accidental" as soon as the symmetry is broken. Un- 
der continuous variations of the circuit's characteristics 
(capacitances and Josephson energies) they move contin- 
uously in the plane (^x = conserving their signature 
and the relation [H, /C] = (but loosing the regularity of 
the conical intersection over the (X, y)-plane). 



IV. MERGING DIABOLICAL POINTS OF A 
COOPER PAIRS PUMP 

As a first approximation in the close vicinity of T^, 
the whole Hilbert space can be reduced to the ground 
eigenspace of HciT"^). In these three- level models, the 
positions of are easily found. They are located in the 
plane ip^ — n at 



,(2?±)«ng(T±)± 



^jr(^jl 



El) 



Suppose that the CPP is mirror symmetric and Ej tun- 
able. The above formula illustrates that the displacement 
of conserves the symmetry [H^V] = 0. Furthermore, 
reducing Ej improves ngi(D~^) = ng^{D^) as much as 
it reduces ngi{D~) = ng2{D~). So the two diabolical 
points _D+ and shown in Fig. 5(b) are expected to 



FIG. 6. The new circuit. The central junction of the CPP 
is replaced by a dc SQUID threaded by a controllable flux 
<l?x. It behaves essentially as an effective junction having a 
capacitance C and a Josephson energy tunable between and 
2 -E.jE . 



merge symmetrically at their midpoint / located at (i, i) 
in the rig-planc. 

Let £^jE = i?,iL = -EjR be the Josephson energy of the 
exterior junctions. By the definition of the capacitance 
unit, their capacitance is ^-j^- Around the point /, an 
approximate Hamiltonian is the restriction of H to the 
subspace spanned by the basis {|0, 0), |1, 0), |0, 1), |1, 1)}. 
Up to an unimportant shift of the zero of energies, the 
truncated Hamiltonian has the form (1) in this basis, 
where X and Y are redefined as follows: 



X 



ng2 + "gi - 1 
1-C 



Y 



1 + 3C 



c 



(1-C)(1+3C) 



= and 



The other parameters are ^ 
G = In this context, hl{Z) and T are the respective 
restrictions of V{ip^)P and e-*('^i+^=)Z^(Z)5'P. Within 
the four-level approximation, if Ej is greater than the 
critical value 



E, 



E^E 

2e 



+ 0{EjE), (7) 



the degeneracies are located at = =F(^d, 0, 0), with 



= - ^) {Ej + AOiEj + Ej, + 40 . 

The merging of the diabolical points and is 
possible if Ej is adjustable. It is well-known that a tun- 
able effective Josephson coupling can be realized via two 
junctions in a loop configuration (a dc SQUID). Such a 
circuit element is de facto interesting to tune couplings 
between superconducting qubits"**^'"*^. It has also demon- 
strated its utility for Cooper pairs pumping in the so- 
called Cooper pairs sluice^®. As shown in Fig. 6, we re- 
place the central junction by a dc SQUID and assume 
all the junctions of the circuit identical. To be consis- 
tent with our previous notations, we set Eje the Joseph- 
son energies and § the capacitances of all the junctions. 
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FIG. 7. Plots of the critical value Ejc as a function of Ej-e, 
obtained by the four-level model on one hand and by a nu- 
merical treatment on the other hand. The plot unit is Ec- 

The new central element has a capacitance C and an 
effective Josephson energy Ej = 2 Eje \ cos ^ | where 
ip'^ = Stt^^. The charging and Josephson Hamiltonians 

still read (3) and (4) after the replacements <^x — > <y5x + ^ 
and U{(px) U{(px,(p'^), the exact definition of the last 
unitary operator being irrelevant for our purpose. We 
also have Cs = 2 C and ^ = 0.4. 

The merging is done by tuning the central coupling 



(through ip[J while we use to maintain the new half- 

fiuxoid condition (fx + ^ = tt. A numerical simulation of 
the process was made, using the 62 closest states of the 
point / to define the truncated Hilbert space. The results 
are in good accordance with the four-level model in the 
Coulomb blockade regime. For example, it is shown in 
Fig. 7 the plot of the critical value Ejc as a function of 
Ej: the numerical result coincides with the expression (7) 
in the limit Ej <^1. 

V. CONCLUSION 

We have demonstrated the possibility of merging dia- 
bolical points of a superconducting quantum circuit. We 
have emphasized the role played by the symmetry for 
that phenomenon. In an experimental perspective, the 
principal difficulty to overcome is the mirror symmetry 
since it is impossible to fabricate two identical junctions. 
This problem can be partially eliminated by using bal- 
anced SQUIDs^^. The theoretical study was accurately 
based on a four-level model whose eigenproblem is ex- 
actly solvable within the constraint of a degeneracy. For 
subsequent works, it may serve as a formal model to im- 
plement mergings of diabolical points in different con- 
texts, such as quantum circuits or cold atoms. It may 
also be used to produce non-Abelian holonomies as well 
as non-Abelian pumpings. 
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